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METHODE POUR DETERMINER DES ZONES DTO MILIEU STRATIFIE, OU 
L'INTERFACE ENTRE UN FLUIDE EN PLACE DANS LE MILIEU ET UN 
FLUIDE DE BALAYAGE, SE DEPLACE DE FAgON STATIONNAIRE 

10 

La pr6sente invention conceme une m6thode pour determiner des zones, d'un milieu poreux 
stratifi6 dans lesquelles le front ou interface entre des fluides en place et des fluides de 
balayage se deplace de fagon stationnaire, c*est-k-dire sans deformation et k vitesse 
constante. 

15 La m6th6de trouve des applications notamment quand on a k balayer par des. fluides 
inject6s (de Teau par exemple) un gisement dliydrocarbures stratifi6. Cette determination 
permet d'obtenir tres rapidement les zones pouvant etre agregees et constituer des unit6s 
hydrodynamiques lors de la mise k r^chelle d*un modele de reservoir. 

Etat de la technique 

20 La simulation des 6coulements diphasiques en milieu poreux heterogene peut requerir 
d'importantes ressoiirces informatiques, en particulier lorsque le modele num6rique du 
milieu consid6r6 est fortement d6taill6. C'est le cas notamment en ingenierie de reservoirs, 
dans le domaine petrolien Afin de pouvoir proc6der k la simulation d'dcoulements avec 
des moyens raisonnables, il faut utiliser une description reduite du reservoir: Une technique 

25 consiste k agreger les mailles du module numerique de fagon k obtenir un maillage plus 
grossier, compose de mailles aux valeurs effectives moyennes. Cette technique, dite des 
pseudos-fonctions, a initialement ete proposee pour des reservoirs stratifies de fa^on k 
agreger les mailles dans la direction verticale par 

- Coats, K.H., Nielsen, RL. and Terhune, M.H. ''Simulation of three-dimensional, two- 
30 phase flow in oil and gas reservoirs.'' SPE 1961, 1967 
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Cette mSthode, reprise et fitendue par de nombreux auteurs, repose sur des hypotheses 
d'^quilibre . vertical capillaire ou d'6quilibre gravitaire oil les forces capillaires ou 
gravitaires sont prddominarites devant les effets dus k la viscosit6, et la distribution des 
saturations dans le reservoir pent etre connue sans recourir k une simulation sur le maillage 
precis (simulation fine). On pent en tirer facilement des valeure moyennes suivant la 
direction verticale. Une seule couche aux valeurs moyennes pent finaleanent suffire k 
d6orire r6volution des saturations dans le r6servoir. 

Une autre m€thode conduisant k obtenir des valeurs moyennes est propos6e par : 

- Heam, C.L. "Simulation of stratified waterflooding by pseudo-relative permeability 
curve. " Journal of Petroleum Technology, pages 805-813, July 1971. 

EUe repose sur I'hypoth^se d'6quilibre vertical dtt aux effets visqueux. Dans ce cas, les 
r^^uilibrages visqueux dans le reservoirs sont tr^s rapides par rapport aux variations de 
saturation, de telle sorte que des valeurs moyennes peuvent toujours 6tre calcul6es sans 
recourir k la simulation fine, car les fluides se propagent k Vitesse constante dans chaque 
couche. La m6thode de Heam, ult^rieurement 6tendue par, entre autres par : 

- Simon, A.D. and Koederitz, L.F. "An improved method for the determination of 
pseudo-relative permeability data for stratified systems." SPE 10975, 1982, 

se base sur un calcul it^ratif conduisant k ordonner les couches du reservoir en fonction de 
la Vitesse de propagation des fluides dans celles-ci. On pent montrer que ce calcul it6ratif 
peut ne pas converger. On est alors en presence d'un fix)nt stationnaiie dans au moins deux 
couches du r&ervoir. Les auteurs proposent alors de ne prendre qu'une seule' couche aux 
propridt6s moyennes pour ces zones particuli^res. Une 6tude th^orique de r6quilibre 
vertical et de ses implications peut 6tte trouv€e chez : 

- Yortsos, Y.C. "Analytical studies for processes at vertical equilibrium." SPE 26022, 
1992. 

Lorsque la capillarity ou la gravity ne peuvent 6tre consid6r6es comme pr^domdnantes sur 
les effets visqueux, mais que l'6quilibre vertical visqueux ne peut se mettre en place, les 
mdthodes propos6es jusqu'a pr&ent ne peuvent se passer d'une simulation fine sur tout ou 
partie du reservoir pour calculer les propri6t6s moyennes pouvant 8tre affect^es aux blocs 
grossiers. On parle alors de pseudo fonctions dynamiques. Celles-ci furent introduites 
notamment par : 
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- ' Kyte, J.R. and Beiry, D.W. "New pseudo functions to control numerical dispersion." 

SPE 5105, 1975. 

Ces mdthodes dynamiques posent toutefois de nomhreux piobldmes th^oriques et 
pratiques. Des reyues de ces mdthodes et des probldmes assocife peuvent 8tre trouvees 
dans les publications suivantes : 

- Archer, R. "Pseudo junction generation". Master of Science thesis, Departement of 
Petiroleum Engmeering of Stanford Univeisity, 1996. 

- Ahmadi, A. "Utilisation des propriStis iquivalentes dans les modules de riservoir : 
cos des Scoulements diphasiques incompressibles." Thdse de Doctoral, Universite de 
Bordeaux 1, 1992 ; 

- Barker J.W. and Thibeau S. "A critical review of the use of pseudo relative 
permeabilities for upscaling". SPE 35491, 1996. 

Les m^tiiodes dynamiques ont en general recours k une agregation systematique et 
uniforme des mailles. Afin de prendre en compte les influences locales des hdt6rog6n6it^s 
sur I'^coulement, des m6tiiodes d'agrdgation non-unifoimes ont 6t6 proposes notamment 
par : 

■ Dannan. N.H. and Durlofsky. L.J. "Upscaling immiscible gas displacements: 
Quantitative Use of Fine Grid Flow Data in Grid Coarsening Schemes." SPE 59452, 
2000. ' 

Ces m6thodes ont Tavantage de n'agr€g©r que des zones pr^f^ientielles et laissent en 6tat 
initial les parties du reservoir k caractferes hydrodynamiques plus complexes. Elles ne 
peuvent cependant s'affranchir d'une simulation fine. 

Des m6thodes permettant d'agr6ger de fagon non unifoimes les strates d'un reservoir, sans 
avoir lecours k une simulation d6taill6e, ont €\6 proposdes. On pent citer des m^tiiodes 
entidrement statiques, comme : 

^- li, D. and Beckner. B. "Optima uplayering for scaleup of multimUlion-cell geologic 
models", SPE 62927, 2000 

ou bas&s sur les r6sultats d'une simulation d'dcoulement monophasique (qui est beaucoup 
rapide qu'flhe simulation d'&oulement diphasique), comme : 
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- Stem, p. and Dawson, A.G.- "A technique for generatirig reservoir simulation grids to 
preserve geologic heterogeneity", SPE 51942, 1999. 

De telles m^thodes, par definition, ne prennent done pas en compte les effets dynamiques 
et.visqueux des ecoulements diphasiques. 

L'6voluUon du fi:ont dans le reservoir au cours de r^coulement est consid6rablement 
influenefe par le couplage visqueux entre le champ de pression et le champ de saturation. 
En particuKer, lorsque le fluide inject^ est moins visqueux et par consequent plus mobile 
au niveau du front que le fluide en place, les instabiUt^ visqueuses favoriseront toujours 
recoulement des fluides dans les couches les plus pemieables. Le temps de percde au 
travers de Celles-ci est beaucoup plus rapide que dans le reste du reservoir. Au contraire, si 
le fluide inject^ est moins mobile, le couplage visqueux peut le ralentir dans les couches 
initialement plus rapides, compensant ainsi les differences de permeability dues k la 
stratification. Un firont stationnaire apparait alors. 

Si, dans au moins une partie du reservoir, il existe un tel front stationnaire, tout se passe 
comme si I'ecoulement n'y voyait qu'une seule couche au sens hydrodynamique. On peut 
done raisonnablement penser qu'une seule couche aux proprietes moyennes dans cette zone 
suffit k la modelisation lors de la simulation numerique de I'ecoulement. Ced est 
d'importance capitale lorsque Ton passe du module geologique de reservoir avec toutes les 
couches, au module de simulation, pour lequel seules les zones les plus influentes sur 
I'ecoulement suffisent k la description. 

Deux probBmes se posent alors. Connaissant la viscosite du fluide injecte, il faut pouvoir 
determiner a priori I'existence et la locaUsation de zones oil I'ecoulement est stationnaire, 
sans recourir pour autant h. une simulation complete d'ecoulement sur le module detailie. D 
faut egalement 6tre en mesure avec la mSme contrainte, de choisir la viscosite du fluide 
injecte de telle sorte qu'un maximum de couches aient un comportement stationnaire. 

La methode que nous proposons permet de reconnaltre les zones d?un reservoir stratifie se 
comportant comme une seule couche au sens hydrodynamique sans recourir k une 
simulation fine et saiis invoquer les hypotheses d'equilibie vertical. Elle faciUte le choix 
des zones k agreger dans les reservoirs stratifies. BUe permet done de prendre en compte 
les effets dynamiques et visqueux tout en permettant une. determination des couches 
grossieres tr&s rapide par rapport aux aneiennes solutions. 
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La m€ithode selon rinvention 

La mdthode selon Tinvention peimet de determiner dans un milieu stratifie dont les 
propri6tes physiques sont connues ou estim^es, au moins une zone oii ime interface entre 
un fluide en place dans le miHeu et un fluide de balayage, de viscosites et densit€s 
diff6ientes connues. injectds dans le miHeu, se deplace de fafon stationnaire, dans le but de 
simpUfier la construction d'un modele de simulation des dcoulements dans le miUeu. Elle 
comporte essentiellement les 6tapes suivantes : 

- on considfere une forme d'interface a priori en supposant que les d6placements de 
fluides en tous ses points sont statiotmaires ; 

- on determine le champ de pression de part et d'autre de I'interface a priori ; 

- on modifieiterativement la forme derinterfacejusqu'kobtenir que les pressions depart 
et d'autre d'au.moinS une partie de I'interface s'dgalisent en tous points de cette partie ; 
et ■ . . ' . 

- on attribue des propri6t& hydrodynamiques moyennes uniform6ment k cha^ue zone du 
milieu delimits par chaque partie d'interface, quand la dite ^galisation est atteinte. 

Suivant un mode de mise en ceuvre, quand on ne pent obtenir une 6gaUsation des pressions 
de part et d'autre de I'interface tout le long de ceUe-ci, on segmente I'interface en plusieurs 
parties et on modifie it6rativement et s6parement la forme de ces diff6rentes parties, jusqu'k 
obtenir une 6galisation des pressions de part et d'autre de chacune d'elles, I'etendue de 
chaque partie d'interface, quand la dite "6galisation est atteinte, d61imitant une zone 
favorable k laquelle on attribue uniform6ment des pK)pri6tes hydrxxiynamiques moyennes. 

Suivant un exemple de mise en ceuvre, on determine la forme d'au moins une zone du 
milieu delimits par une interface k dgplacement stationnaire, qui correspond k des valeurs 
differentes de la viscositd des fluides de balayage. et on choisit la viscosity pour laquelle on 
optimise les d6placements stationnaires dans le dit milieu. 

Pr&entation succincte des figures 

Les caracteristiques et avantages de la m6thode selon invention, apparaittont plus 
clairement k lar lecture de la description ci-apr^ d'un exemple non limitatif de mise en 
ceuvre. en se r6f(6rant aux dessins annexes oCl : 
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- la figure 1 montre une boucle de d^tennination it^rative de I'existence et de la fonne 
d'un front stationnaire pour un 6coulement diphasique de rapport de viscosit6s connu 
dans un milieu stratifie donn^ ; 

- la figure 2 montre une boucle de construction it^tive de revolution des couches 
hydrodynamiques en fonction du rapport des viscosites, pour des 6coulements 
diphasiques dans un milieu stratifi6 donn6 ; 

- la figure 3. montre un exemple de r6servoir stratifi6 incline par rapport k la verticale. 
Cette figure a pour but d'expliciter les notations employees dans la description de la 
m6thode; 

- la figure 4 montre la valeur des perm6abmt6s absolucs en fonction de la profondeur. 
employees dans rexenq)le synth6tique de i^ervoir illustiant la m6thode de fagon non- 
limitative ; 

- la figure 5 montre revolution de I'interface au cours du processus it&atif de recherche 
de la forme d'un fix>nt stationnaire, sur I'exemple d6ciit par la figure 4 ; 

- la figure 6 montre revolution de la solution stationnaire en fonction du rapport des 
viscosites employees, sur I'exemple deciit par la figure 4 ; ct " 

- la figure 7 montre revolution du nombre de couches hydrodynamiques en fonction du 
• rapport des viscosites employees, sur I'exemple decrit par la figure 4. 

Description d^taill^ de la methode 

On considere un balayage de Guides de viscosites connues dans un mlKeu stratifie. 
Supposons qu'un front stationnaire de forme connue se developpe dans le miHeu consideie. 
Ceci impose une condition de flux au niveau de I'interface. Le flux est egalement connu sur 
les bords du reservoir '(oii U est en general nul), en entree et en sortie (oh il est impose par 
les conditions d'injection et de lecupetation). Comme la pression est une variable continue, 
le saut de pression de part et d'autre de I'inteiface est nul en tdut point de celle-ci. 

Cta cherche k present I'existence et la forme d'un front stationnaire dans le milieu. On 
procede de fa^on iterative en partant d'une forme de front a priori, en supposant que celle- 
ci correspond k la solution stationnaire. Comme le front est suppose statiomiaire, on 
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' connait le flux k Hnterface et Ton a done un probBme de Neumann parfaitement d6fini k 
i^soudie pour chacune des deux zones que separe I'interface. Le champ de pression peut 
done etre r6solu ind6pendaniment dans les deux regions. On lvalue alors le saut de 
pression de part et d'iautre de I'interface, en tout point de ceUe-ci. Si le saut de pression est 
nul en tout point, on possede alors la forme correcte du front coirespondant k un 
6coulement stationnaiie. Si le saut n'est pas nul partout, la condition de continuity des 
pressions n'est pas respect6e et la forme de front ne correspond pas k une solution 
stationnaire. On doit alors d6fonner I'interface et recommencer le processus jusqu'k la 
convergence. Deux cas apparaissent alors : 

- Si la convergence est atteinte, on possMe alors la forme du front conespondant k une 
solution stationnaire de I'ecoulement. Toutes les couches sont traverses par le fluide k 
la meme vitesse, et I'ensemble du systeme se comporte comme une seule unit^ 
hydrodynamique. On peut done ne consid6rer qu'une seule couche aux propri6t6s 
moyennes lors de la mise k I'dchelle. 

- Si le processus it^ratif ne converge pas, e'est qu'un front stationnaire ne peut pas exister 
sur I'ensemble du milieu. Lors de la deformation du front, il apparait alors en g6n&al 
deux zones s'61oignant perp^tuellement l-une de I'autre. On d6coupe alors le milieu en 
deux parties correspondants aux zones pr6c6denQment identifies. Puis le processus de 
recherche d'une solution stationnaire reprend dans chacune de ces deux parties. En 
proc&iant ainsi de mani^re iterative, on finira par obtenir toutes les zones dans 
lesquelles des fronts stationnaires peuvent se d^velopper, et qui peuvent Stre 
consid6r6es comme des unites hydrodynamiques. 

La m^thode peut aussi 6tre appUqu6e de fagon iterative en faisant varier le rapport des 
viscosity des fluides consid^r^s. On part dhm fluide inject6 beaucoup plus visqueux que 
celui en place. On doit alors trouver rapidement une solution stationnaire sur toute 
l'6paisseur du miUeu. Lorsque Ton diminue le report des viscosity entre le fluide mject6 
et le fluide en place, la forme du front stationnaiie, si celui-ci existe, doit 8tre modifife. On 
reprend alors la mdthode ddcrite ci-dessus pour trouver la nouvelle forme du front, en 
utilisant comme forme de front initiale du processus itdratif, la forme finale obtenue pour 
le precedent rapport des mobility. Ceci permet de converger beaucoup plus rapidement 
vers la solution stationnaire. 
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On peut ainsi procMer de fa9on iterative sur le rapport des viscosit6s. Au fur et ^ mesure 
que Ton d&rolt ce rapport, le couplage visqueux ne sera plus suffisant pour compenser les 
contrastes de penn6abiUte et une solution stationnaire ne pouira plus 6tre trouvfe sur toute 
r^paisseur du miUeu. On d6coupe alors le milieux en deux portions suivant la mfithode 
d^crite ci-dessus en cas de non-convergence, et la recherche se pouisuit dans les deux 
portions de miUeu. En proc^t ainsi de fagon it^ve sur le rapport des viscosit^s, on 
obtient le nombre de couches hydrodynamiques du noiUeu en fonction des viscosit^s 
considdr^es. Ce i^sultat pourrait pennettre par exemple de connattr^ la viscosity du fluide 
d'mjection qui pMmet d'optimiser la r6cup6ration. 

On considdre une section verticale dHm reservoir stratifie selon une direction x inclin^e par 
rapport k ITioiizontale (figure 3). valeurs des penn6abiUt6s absolues sont connues dans 
toutes les strates. L'6paisseur du reservoir est selon la direction z. A I'extrgmit^" amont {x = 
0) on injecte. k d^bit constant Qo impos6, un fluide incompressible de viscosity m et de 
density qui pousse de fa9on immiscible le fluide initialement present dans le reservoir, 
incompressible et de Viscosity ^/z et .de densit6 p2. On suppose qu'il n'y a pas de pression 
capiUaire entte les deux fluides. La recuperation se fait a rextremitd aval {x = L). Le flux 
est nul i travels le toit = 0) et le mur du reservoir {z = A). On suppose egalement que les 
courbes de penn^abilites relatives sont uniformes dans le reservoir et que chacune de ces 
deux zones est k saturation constante en fluides considei^s. La zone non encore balay^e est 
done h saturation irr^ductible Si, en fluide inject^ et la zone balayde k saturation maximuni 
S„ on saturation de front Sf selon les hypotheses de IHitiHsateur. Le front de saturation est 
arbitrairement plac6 au centre du rgservoir, avec une forme pr^supposfe qui correspond k 
notre a priori sur la forme du firont stationnaire. Cette forme sera par la suite modifiee 
jusqu'k trouver la "vraie" solution stationnaire, si ceUe-ci existe. On doit done r^soudre un 
probldme aux firontieres libres. 

Les conditions aux limites et la condition de front stationnaire pennettent de definir un 
probieme de Neumann dans les deux zones que sSpare le front : 

- A I'extrimiti amont : la condition d'injection impose un flux constant sur I'epaisseur du 
reservoir en entree. Si le milieu est suppose suffisamment long et le profil de 
saturations deji bien developpe. on peut supposer que le gradient de pression en entree 
du systeme est identique.dans toutes les couches : 
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dx 



x=0 

dp 

dz 



- >"igo 



2=0 



ohp est la pression au point consid6r6 et est la moyenne des perm6abiUt6s effectives en 
entree siir toute I'^paisseur du milieu, et j2o est le da)it d'injection. Les penn6abiUt6s 
effectives prennent en compte I'effet des perm6abiUt^s relatives, en fonction des saturations 
choisies par I'utilisateur dans les zones amont et aval. 

- A Vextrimiti aval : la condition d'incompressibiUt6 impose que le flux sortant du 
systfeme soit le mgme que le flux entrant! On obtient done une condition simUaire k la 
tranche de r6cup6ration : 



dp 

dx 



Bz 



M2Q0 



oH est la moyenne des peim6abilit6s effectives en sortie sur toute I'epaisseur du mitieu. 



de part de d'autre de I'interface : Comme le front est suppos6 stationnaire. il se d6pli 
horizontalement sans d^fonnation k la vitesse de filtration , d^finie telle que : 



ace 



u. 

A 

Oil X est le vecteur unitaire dans la direction horizontale. Cette condition impose le flux et 
done le gradient de pression de part et d'autre de Tinterface : 



.n 



oii nest le vecteur unitaire orthogonal h I'interface au point consid6r6, Vp, est le gradient 
de pression k I'interface du c6t6 amont et Vp^ est le gradient de pression a interface du 
c6td aval. 
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- Equation en pression : X'6quation en pression est donn6e classiquement -de part et 
d'autre de I'interface. Pour la zone aval : 



Pour la zone amont : 



i 
i 



= 0 



— (P2g-Vp) 

Les ^uations en pression et les conditions aux limites pennettent de d^finir pour chaque 
zone de part et d'autre du front un probl^me de Neumann. Le champ de pression dans 
chaque zone peut done 6tre calculd num&iquement ^ une constante pi^s pour chaque zone. 
On ^limine ce probldme en posant la condition supplementaire suivante, peimettant 
d'61iminer les constantes en recalant les pressions le long du front : 

fp,(z)dz = jpMdz = C 

oiipi et/72 sont les pressions 6valu6e au front respectivement pour les zones amont et aval, 
et C est d6finie arbitrairement par Ihitilisateur et peut 6tte nulle. 

Une fois la pression p calcul6e dans chaque region et recalee le long de I'interface, le saut 
de pression (z) h trav^s le front en tout point de ce dernier peut 8tre calcul6 : 

(z) = P2 (z) r-Pi (z) - A/>(g.x + gje) 

Si le saut de pression est partout nul le long du front, ou du moins en degk de la valeur 
critique d6finie par I'utilisateur, on dispose de la forme du front correspondant k un 
ddplacement stationnaire. Sinon, il faut d^former I'interface de fagon k faire diminuer les 
sauts de pression jusqu'k obtenir la forme stationnaire. Ceci revient alors k r6soudre un 
probldme d'optimisation, pour lequel de nombreuses m^thodes sont applicables. 

Ce processus de minimisation peut ne pas converger. Dans ce cas, il n'existe pas de front 
stationnaire sur toute I'dpaisseur du milieu. La m6thode ddcrite doit done Stre reprise sur 
des parties plus petites du reservoir. On peut ainsi trouver des fronts stationnaires dans 
plusieurs zones du reservoir, qui se d^placent k des vitesses diff6rentes. 
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Exemple ■ , 

On d^t ici un exemple synth6tique iUustrant la m6thode de fagon non-limitative. 

On consid^re un reservoir stratifi6 avec quatre strates hoiizontales de peim6abilit6s 
d.ff6rentes. indiqu^es sur la figure 4. Pour simplifier au maximum, on h6gligera les effets 
de la gravite et on consid^ que les penn6abilit6s relatives des strates sont telles que la 
zone non balayee est satur^ en huile seule et qu'en amont de interface, seule I'eau est 
pr^sente k saturation 5=i. La porosity est constante dans le reservoir. Ueau est inject^e sur 
la face amont h. flux constant connu. LTiuile est h. viscosity de 1. 

On commence par chercher I'existence et la forme d'un front statiomiaiie dans le miUeu 
lorsque la viscosity de I'eau est de 12. Comme on ne dispose pas encore d'a priori, on part 
d'une interface plane, que Ton d6forme de fa^on it6rative jusqu'k atteindre la convergence 
vers une solution stationnaire (figure 5). 

Une fois la solution obtenue pour cette viscosity, on pent rechercher la solution stationnaire 
pour une viscosity de I'eau inf^ieure, par exemple une viscosity de 10. La solution 
statiomxaire obtenue pr&6demment peu servir d'a priori en enlx^ du processus it^ratif de 
d6fonnation du front. On converge alors vers mie nouvelle solution. En proc6dant ainsi 
il^tivement sur diff^ntes valeur. de la viscosity de I'eau, on obtient les formes du front 
statiom^aire associ^es k ces valeurs (figure 6). Comme le rapport des mobilit6s entre le 
flmde ne place et le fluide inject6 diminue chaque fois que I'on diminue la viscosite de 
reau. le front s'^tale au fur et k mesure, car les effets visqueux ont de plus en plus de mal k 
compenser les differences de permeability d'une strate k I'autre. Lorsque la viscosity de 
reau diminue en dega de la valeur critique 5. on s'apergoit que le processus de recherche 
d'un front stationnaire ne converge plus. Les positions du front dans les deuxi^me et 
troisi&me strates s'61oignent perpetueUement 1-un de I'autre. On consid6re alors qu'un front 
stauonnaire ne pent se former sur toute I'^paisseur du milieu. Ce dernier doit 6tre d6coup6 
en deux parties, la separation se situant entxe la deuxifeme et la troisifeme couche On 
poursuit ensuite le processus de recherche de fr^nt statiomiaire dans chacune des parties 
pr&6deinment identifies. 

En iterant ainsi sur la valeur de la viscosite du fluide injecte, on pent construire un " 
diagramme representant la r^ponse du milieu du point de vue hydrodynamique en fonction 
du rapport des mobilites M entre le fluide en place et le fluide injecte (figure 7) On voit 
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sur ce diagramme que lorsque le rapport des mobilit^s est sup6rieur k 5, le milieu se 
c6iiq)orte cdmme une seule couche hydrodynamique, puisqu'un front stationnaire existe sur 
toute son 6paisseur. Lorsque M passe entre 4,5 et 5 deux couches hydrodynamiques 
apparaissent : la premiere est fonn6e des deux premieres couches "g6ologiques", et la 
seconde des deux autres couches "geologiques". Puis les deux premieres couches 
geologiques se s6parent au sens hydrodynamique en degh d'un rapport M de 4,5. On a alois 
trois couches au sens hydrodynamique dans le milieu. Les deux demidres couches 
geologiques se sdpaneait en 6eqk de M=3. On a alors quatre couches hydrodynamiques, qui 
correspondent aux quatre couches geologiques pr&entes dans le miheu. 
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-REVENDIGATIONS 

1) M6thode pour determiner dans un miUeu stratifi6 dont Jes propri6tds physiques sont 
connues ou estimees, au moins une zone oil une interface entre un fluide en place dans le 
milieu et un fluide de balayage, de viscosit^s et densit6s diff&entes connues, inject^s dans 
le miUeu, se ddplace de fa9on stationnaire, dans le but de simplifier la construction d'un 
module de simulation des feoulements dans le miKeu, caract6iis6e en,ce que : 

- on considfere une forme d'interface a priori en supposant que les d6placements de 
fluides en tous ses points sont stationnaiies ; 

- on d6tenmne le champ de pression de part et d'autre de I'interf ace a priori ; et 

- on modifie it^rativement la forme de I'interf ace jusqu'^ obtenir que les pressions de part 
et d'autre d'au moins une partie de I'interface s'egalisent en tous points de cette partie ; 
et 

- on attribue des propri6t6s hydrodynamiques moyennes uniform6ment k chaque zone du 
miUeu d61imit6e par chaque partie dlnterface, quand la dite 6gaUsation est atteinte. 

2) M6thode selon la revendication 1, caractdrisde en ce que, k d6faut d'obtenir une 
6galisation des pressions de part et d'autre de I'interface tout le long de c?elle-ci, on 
segmente I'interface en plusieurs parties et on modifie it6rativement et s6par6ment la forme 
de ces diff^tes parties, jusqu'k obtenir une 6galisation des pressions de part et d'autre de 
chacune d'elles, I'^tendue de chaque partie d'interface, quand la dite dgalisation est atteinte, 
dfilimitant une zone favorable k laqueUe on attribue uniform^ment des propri6t6s 
hydrodynamiques moyennes. 

3) Mgthode selon la revendication 1 ou 2, caract^sfe en ce que I'on d6teiimne la forme 
d'au moins une zone du miKeu d^Umitfe par une interface k d6placement stationnaire, qui 
correspond k des valeurs diff6rentes de la viscosity des fluides de balayage, et on choisit la 
viscosity pour laquelle on optimise les d6placements stationnaires dans le dit miUeu. 



wo 2004/068176 



PCT/FR2004/000077 



1/2 

FIG.1 



X 



1 



FIG.2 



it, ^ 




^6a <Uaiv'(agerice 



accept^. conuricf 



Constiucfion des coudies 
hydro^fnamtqiies en fbncli cm 



wo 2004/068176 



PCT/FR2004/000077 



2/2 



FIG.3 




FIG.7 




1 .2 3 4 5 '6 7 S 9. 



INTERNATIONAL SEARCH REPORT 



atlonal Application No 

/FR2004/000077 



4v9!-ASSlRCATION OF SUBJECT MATTER 

IPC 7 601V11/00 



According !o Internatronal Patent Classlficaiion (IPC) or to both national classification and IPC 



B. RELDS SEARCHED 



Minimum documentation searched (dassiflcallon system foBowed by classification symbols) 

IPC 7 GOIV 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched" 



Eteclronic data base consulted during the international search (name of data base and. where practical, search terms used) 

EPO-Internal , WPI Data 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Categoiy ' 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



SHELDON ET AL.: "The Motion of an 
Interface Between Two Fluids In a slightly 
dipping Porous Medium" 

SPE, SOCIETY OF PETROLEUM ENGINEERS INC., 

1962, pages 275-282, XP002261168 
page 275 -page 280 

GLIMM J. ET AL.: "Front Tracking for 

Petroleum Reservoir Simulation" 

SPE, SOCIETY OF PETROLEUM ENGINEERS INC., 

1983, pages 41-49, XP002261169 
the whole document 



-/- 



1-3 



1-3 



I X| Further documsnis are listed In the conHnuaUon of box C. 



□ 



Patent family members are listed in annex. 



" Special categories of cited documents : 

•A' document defining the general state of the art which Is not 
considered to be of particular relevance 

■E' earlier document but published on or after the international 

filing date 

■L* document which may throw doubts on priority clalm(s)or 
which is cited to estabHsh the pubHcatton date of another 
citation or other special reason (as specified) 

'O" document referring to an oral disdosuiB, use. exhibition or 
other means 

■P' document published prior to the inlernatlonai fiBng date but 
later than the priority date claimed 



Date of the actual completion of the International search 



16 June 2004 



T' later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underiyinq the 
invention 

•X' document of particular relevance; the claimed invention 
cannot be considered novel or cannot l>e considered to 
Involve an inventive step when the document is tai<en alone 

•Y' document of particular relevance; the claimed invention 

cannot be considered to involve an Inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 

document member of the same patent family 



Date of mailing of the International search report 

07/07/2004 



Name and mailing address of the ISA 

European Patent Office. P.B. 5818 Patenllaan 2 
NL-2280HVRljswiik 
TeL (+31-70) 340-2040, Tx. 31 651 epo nl. 
Fax: (+31-70) 340-3016 



Authorized officer 



Schnelderbauer, K 



f=bmi PCT/ISA/210 (socond sheet) (January 2004) 



INTERNATIONAL SEARCH REPORT 



li^Tiatlonal Application No 

Wr/FR2004/000077 



C.(Contlnuation) DOCUMENTS CONSIDERED TO BE RELEVANT 



Category * Citation of document, with Indication, where appropriate, of the relevant passages 



Relevant to claim No. 



YORTSOS: 
Processes 
SPE, SOCIETY 
1992, pages 
cited in the 
page 1 -page 



"Analytical Studies for 
at Vertical Equilibrium" 
OF PETROLEUM ENGINEERS 
1-14, XP002261170. 
application 
12 



1-3 



INC. 



DARMAN ET AL.: "Upscaling Immiscible Gas 
Displacements: Quantitative Use of Fine 
Grid Flow Data in Grid Coarsening Schemes" 
SPE, SOCIETY OF PETROLEUM ENGINEERS INC., 
2000, pages 1-14, XP002261171 
cited In the application 
page 1 -page 14 



1-3 



Fbim PCT/lSA/210 (continuation of second sheet) (January 2004) 



RAPPORT DE RECHERCHE INTERNATIONALE 



Iande Internatlonaie No 

/FR2004/000077 



A. CLASSEMENT DE L'OBJET DE LA DEMANDE 
CIB 7 GOlVll/00 



Seton la dassificalton Internationale des brevets (CIB) ou & la fois selon !a classification nationale et ia CIB 



B. DOMAINES SUR LESQUELS LA RECHERCHE A PORTE 



Documentation mmlmale consuIt6e (sysl6me de classification suivi des symbotes de dassement) 

CIB 7 GOIV 



Documentation consull^e autre que la documentaUon mInlmaJe dans la mesure oCi ces documents relfevent des domaines sur lesquels a port6 la recherche 



Base de donn6es 6lectronlque consult6e au oours de la recherche Internationale (nom de la base de donn6es. et si realisable, termes de recherche uUllsesr 

EPO-Internal , WPI Data 



C, DOCUMENTS CONSIDERES COMME PERTINENTS 



Cat6gorfe ° Identification des documents cfl6s. avec, le cas 6ch6ant. rindication des passages pertinents 



no. des revendicallons vls^es 



SHELDON ET AL.: "The Motion of an 
Interface Between Two Fluids in a slightly 
dipping Porous Medium" 

SPE, SOCIETY OF PETROLEUM ENGINEERS INC., 

1962, pages 275-282, XP002261168 
page 275 -page 280 

GLIMM J. ET AL.: "Front Tracking for 

Petroleum Reservoir Simulation" 

SPE, SOCIETY OF PETROLEUM ENGINEERS INC., 

1983, pages 41-49, XP002261169 
le document en entler 



1-3 



1-3 



j X| ^^^^ ^ O pour la fin de la Qste des documents 



j j Les documents de families da brevets sont indiqu6s en annexe 



Categories sp^clales de documents clt^s: 

■A" document d6flnlssant l'6tat g6n6ral de la technique, non 
consid6r6 comme partlculiferement pertinent 

document ant6rieur, mals public k la date de d6p6t international 
ou apr&s cette date 

■f document pouvant Jeter un doute surune revendicatlon de 
priority ou clt6 pour d6lemilner la date de publication d'une 
autre citation ou pour une ralson sp6cia!e (telle qu*indlqu6e) 

'O* document se r^f^rant k une divulgation orale, k un usage, k 
une exposition ou tous autres moyens 

•P" document pubI16 avant la dale de d6pdt Inlematlonal, mals 
post6rieurement h la date de priority revendlqu^e 



■T" document ult6rieur publl6 aprfes la date de d6p6t international ou la 
date de priority et n'appartenenant pas k r^at de la 
technique pertinent, mals cltd pour comprendre le princIpe 
ou la th^orie constltuant la base de rinventlon 

■X* document partlculiferement pertinent; rinven tion revendlqu6e ne peut 
£tre consld^r^e comme nouvelie ou comme impliquant une activity 
Inventive par rapport au document consider^ isoldment 

•V document particullferement pertinent; rinven lion Fevendiqu^e 
ne peut 6tre consid6r§e comme bnpDquant une activity inventive 
(orsque le document est assode k un ou plusieurs autres 
documents de m6me nature, cette comblnaison 6tant dvidente 
pour une personne du metier 

document qpl fait partle de la mdme famllte de brevets 



Date a iaquelle la recherche Internationale a M effecUvement achev6e 

16 juin 2004 



Dale d'exp^ditlon du present rapport de recherche Internationale 

07/07/2004 



Nom et adiesse postals de radministrallon charg6e de la recherche Internationale 
Office Europ^en des Brevets, P,B. 5816 Patentlaan 2 
NL-2280HVRIlswlik 
Tel. (+31-70) 340-2040. Tx. 31 651 epo nl. 
Fax: (+31-70) 340-3016 



Fonctionnalre autorisd 



Schneiderbauer, K 



Formulairo PCT/ISA/210 (deuxi&me feuilla) (Janvier 2004) 





RAPPORT DE RECHERCHE INTERNATIONALE 1 

1 Wr/FR2004/000077 




C.(suite) DOCUMENTS CONSIDERESCOMME PERTINENTS 




Catfgorie ' 


Identification des documents ctt^, avec, le cas ech^ant, Tindlcation des passages pertinents 


no. des revendications vis^es 




A 


YORTSOS: "Analytical Studies for 

D«M A MAMA If 4 M 1 ^ ■ ■ M ^ H Lk. (A m ■ II 

processes at vertical hqumDnum 

SPE, SOCIETY OF PETROLEUM ENGINEERS INC., 

1992, pages 1-14, XP002261170 
cite dans la demande 
page 1 -page 12 


1-3 




A 


DARMAN ET AL.: "Upscaling Immiscible Gas 
Displacements: Quantitative Use of Fine 
Grid Flow Data 1n Grid Coarsening Schemes" 
SPE, SOCIETY OF PETROLEUM ENGINEERS INC., 

2000, pages 1-14, XP002261171 
cite dans la demande 
page 1 -page 14 


1-3 



Fonnuialre PCT/ISA/&10 (suite de la demdtoM feuDle) ganvier ao04) 



